Quantum Mechanics of Damped Systems 

Dariusz Chruscinski* 
Department of Mathematics and Statistics 
University of Calgary 
Calgary, Alberta, Canada 



Abstract 

We show that the quantization of a simple damped system leads to a self-adjoint 
Hamiltonian with a family of complex generalized eigenvalues. It turns out that they 
correspond to the poles of energy eigenvectors when continued to the complex energy 
plane. Therefore, the corresponding generalized eigenvectors may be interpreted as res- 
onant states. We show that resonant states are responsible for the irreversible quantum 
dynamics of our simple model. 
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1 Introduction 



Standard textbooks on quantum mechanics investigate mainly the Hamiltonian system, i.e. 
by a quantum system one usually means a Hilbert space TC which describes physical quantum 
states and a self- adjoint operator (Hamiltonian) in TC which governs dynamics of the system. 
However, most of the classical systems are not Hamiltonian and the quantum mechanics of 
such systems is poorly understood. In the present paper we are going to investigate one of 
the simplest non-Hamiltonian system corresponding to a damped motion in one dimension: 

x = — 7X , (1-1) 

where x G R, and 7 > stands for the damping constant. Classically, the damping behavior 
is described by the exponential law 

x(t) = e'^x . (1.2) 

As is well known [1] (see also [2]), within the standard Hilbert space formulation of quantum 
mechanics there is no room for such a behaviour on a quantum level. Therefore, in order to 
deal with this problem, we shall use the rigged Hilbert space approach to quantum mechanics 
which generalizes the standard Hilbert space version [3, 4, 5]. A rigged Hilbert space (or a 
Gelfand triplet) is a collection of spaces [6, 7]: 

$CHC$', (1.3) 

where TC is a Hilbert space, $ is a dense nuclear subspace of TC, and $' denotes its dual, i.e. 
the space of continuous functionals on <E> (see section 2 for a brief review). 

The quantization of our simple model (1.1) leads to a self-adjoint Hamiltonian H in 
TC = L 2 (IR). Interestingly, H being self-adjoint, gives rise to the family of generalized complex 
eigenvalues. Clearly, these eigenvalues are not elements of the spectrum cr(H) = (—00,00). 
The corresponding eigenvectors do not belong to L 2 (R) but to $' for an appropriately chosen 
<!>. We show that these complex eigenvalues have many remarkable properties analogous to 
the point spectrum of a self-adjoint operator. In particular, they give rise to the spectral 
decomposition of H. Moreover, they are closely related to the continuous spectrum of H. It 
turns out that they correspond to the poles of the energy eigenvectors ip E when continued to 
the complex energy plane [8]. Physicists usually called the corresponding eigenvectors reso- 
nant states [3, 9, 10] (see also [11]). It is widely believed that resonant states are responsible 
for the irreversible dynamics of physical systems (see e.g. recent collection of papers [3]). 
Indeed, it is true in our simple model. To see this we construct two Gelfand triples: 

$± C L 2 (R) C $' ± , (1.4) 

such that <J> + n<I>_ = {0}. Obviously, the time evolution is perfectly reversible when considered 
on L 2 (R). It is given by the 1-parameter group of unitary transformations U(t) = e~ tHt . 
However, when restricted to <3?±, it defines only two semigroups: U(t > 0) on <&_, and 
U(t < 0) on $_|_. Therefore, the evolution on $± is irreversible. This irreversibility is caused 
by quantum damping, or, equivalently, by the presence of resonances. 
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2 Rigged Hilbert space 

Consider a rigged Hilbert space, i.e. the following collection (Gelfand triplet): 

$cHcf', (2.1) 

where TC is a Hilbert space with the standard norm topology t-h, is a topological vector 
space with a topology, r$, stronger than t%, and <£>' is the dual space of continuous linear 
functionals on <E> [6, 7]. We denote the action of <£' on <3? using Dirac notation, i.e. for any 
<p € $ and F6$' 

<#F):=F(<£) . (2.2) 
Any self-adjoint operator A in Tl may be extended to an operator on 

(2.3) 

by 

(0|lF):=(^|F). (2.4) 

Now, if for any ^ G $ 

(cf>\AF x ) = \{<f>\F x ) , (2.5) 

then Fx G <&' is called a generalized eigenvector corresponding to a generalized eigenvalue A. 
Omitting (j) one simply writes: 

A\F X ) = \\F X ). (2.6) 

Note, that a generalized eigenvalue A may be complex. Now, if the spectrum of A 

a(A) = a p (A) U a c (A) C R , (2.7) 

with <7p(A) = {Ai,A2,...}, then the Gelfand-Maurin theorem [6, 7] implies the following 
spectral decompositions: 



U = Y,\F n )(F n \+ f ^d\\F x ){F x \, (2.1 



and of A itself: 



A = y^\ n \F n ){F n \+ [ d\\\F x )(F x \. (2.9) 

This way the rigged Hilbert space approach fully justifies the standard Dirac notation. 

The choice of depends on the particular problem one deals with. In the present paper 
we shall consider the following functional spaces: V - the space of C°°(R) functions with 
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compact supports equipped with the convex Schwartz topology [12], S - the space of C°°(R) 
functions vanishing at infinity faster than any polynomial [12]. Moreover, let us define 

Z:={F[^]\^eV} , (2.10) 

where F[(j>] denotes the Fourier transform of (f>. It turns out [13] that Z is isomorphic to the 
space of entire functions of fast decrease along R. More precisely, let 

F L [<f>](z) := -4= f°° e izx cP(x)dx , (2.11) 



J- 

be the Fourier-Laplace transform of <fr G T>. One proves [12, 13] the following 

Theorem 1 (Paley- Wiener- Schwartz) Let a > 0. An entire function U(z) is a Fourier- 
Laplace transform of a function u £ V with support 

supp(u) = {x G R | \x\ < a} , 

if and only if 

\z\ n \U{z)\ < C n e a ^ lmz \ , n = 1,2,... . 
Now, for z = x G R, i.e. Imz = 0, Fl[4>] = F[4>], and the above theorem implies 

\x\ n \F[4>](x)\ <C n , n= 1,2,... . (2.12) 
Clearly, Z f]T> = {0}. Moreover, one has 

VcS C L 2 (R) , (2.13) 

and 

Z cS C L 2 (R) , (2.14) 

and both V and Z are dense in S. One proves [13] that the Fourier transformation which 
defines the unitary operator 

F : L 2 (R) — ► L 2 (R) , (2.15) 
establishes an isomorphism between V and Z. 

3 Quantization of damped systems 

Let us quantize a classical damped system described by (1.1). Clearly this system is not 
Hamiltonian. However, it is well known (cf. [14]) that any dynamical system may be rewritten 
in a Hamiltonian form. Consider a dynamical system on n-dimensional configuration space 
Q: 

x = X(x) , (3.1) 
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where X is a vector field on Q. Now, define the following Hamiltonian on the cotangent 
bundle V = T*Q: 

H(a x ) :=a x (X{x)) , (3.2) 
where a x £ T*Q. Using canonical coordinates (x 1 , . . . , x n ,p\, . . . ,p n ) one obtains: 

n 

H(x,p) = Y,PkX k {x) , (3.3) 

k=l 

where X k are components of X in the coordinate basis d/dx k . The corresponding Hamilton 
equations take the following form: 



x k = {x k ,H}=X k {x) , (3.4) 

)X\x 
dx k 



p k = {Pk,H} = -± Pl ^-, (3.5) 



l=i 

for k = 1, . . . , re. In the above formulae { , } denotes the canonical Poisson bracket on T*Q: 

k=l 

Clearly, the formulae (3.4) reproduce our initial dynamical system (3.1) on Q. 

Let us apply the above procedure to the damped system (1.1). One obtains for the 
Hamiltonian 

H(x,p) = —~fxp , (3-7) 
and hence the corresponding Hamilton equations 

x = -jx , P = 1P , (3.8) 
give rise to the following Hamiltonian flow on M 2 : 

{x,p) — ► (e~ 7 *x,e'» . (3.9) 

Now, the quantization is straightforward: one has for the Hilbert space TL = L 2 (M), and for 
the Hamiltonian 

H = -^(xp + px) . (3.10) 

It is evident that (3.10) defines a symmetric operator on L 2 (IR). In section 4 we show that 
H is self-adjoint and hence it gives rise to a well defined quantum mechanical problem. 1 

1 Actually, this Hamiltonian is well known in quantum optics in connection with the squeezed states of light 
[15]. Introducing a and a*: 

„ a + a* „ a — a* 

x = -=- , P ■- 



v/2 ' 1 V2i ' 
the Hamiltonian (3.10) may be rewritten as follows: 

H= Y .(a -a) , 

which is exactly a generator of squeezing. 
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Let us observe that performing the canonical transformation 



(-yX-P), p=-$={ 1 X + P) , (3.11) 



the classical Hamiltonian (3.10) takes the following form: 

H = \{P 2 ~1 2 X 2 ) , (3.12) 

that is, it corresponds to the so called reversed harmonic oscillator. This system was analyzed 
in [16] and recently in [17, 18, 19] (see also [20, 21]). 

4 Properties of the Hamiltonian 

Let us investigate the basic properties of the Hamiltonian defined in (3.10). 
Proposition 1 The operator H = — ^(xp + px) is self-adjoint in L 2 (M.). 

Proof. To prove that H is self-adjoint we show that e~ lH is unitary in L 2 (M.). One has 

H = --ixp + px) = i7 ( x— + - j . (4.1) 



Let us define 



Clearly, 



U = e~ iS = e ^ /2 e lxdx . (4.2) 



U^(x) = e 7/2 VKe 7 x) , (4.3) 
for any ip G L 2 (R). The operator U defines an isometry: 



(Uijj\U4>) = / Uip(x)U(f>(x)dx = / e 1 ^{enx)4>{e l x)dx= I ^(y)^(y)dy 

J -co J — oo J — oo 

= <V#>- (4.4) 

Moreover, due to (4.3), U is onto, and hence it is unitary in L 2 (R). Therefore, Stone's 

theorem implies that H is self-adjoint (see e.g. [12]). □ 
Obviously, H is parity invariant: 

PHP' 1 = H , (4.5) 
where the parity operator P is defined by: 

PxP" 1 = -x , PpP 1 = -p . (4.6) 
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Now, let us turn to the time reversal operator T. The theory invariant under the time reversal 
has the following property: if ip(t) is a solution of the Schrodinger equation given by 

m = U{t)ij> , (4.7) 

with U(t) = e~ iHt , then Tip evolves into 

(Tty)(-t) = l7(t)(TV0 , (4.8) 

or, equivalently 

T(l7(t)V0 = U(-t)(T1>) , (4.9) 

for any tp £ TC. Now, following Wigner [22], T is either unitary or anti-unitary. If T is 
unitary, then (4.9) implies 

TH + HT = 0. (4.10) 

It means that if 

Hip E = Eip E , (4.11) 

then 

HTip E = —ETip E , (4.12) 

that is, any eigenvector ip E with the energy E is accompanied by Tip E with energy —E. 
Usually, this case is excluded since one expects that the Hamiltonian is bounded from below. 
If this is the case, then T is anti-unitary and (4.9) implies: 

TH-HT = Q. (4.13) 

However, the Hamiltonian defined in (3.10) is not bounded from below, and, as we show in 
section 6 its spectrum a(H) = (-co, oo). Therefore, we take T to be unitary in L 2 (R). 

Proposition 2 The time reversal operator T is realized by the Fourier transformation: 

Tip := F[ip] , (4.14) 

i.e. 

F^HFip = -Hip , (4.15) 

for all ip G L 2 (M). Moreover, 

T 2 ip(x) = P<f>(x) = ip(-x) . (4.16) 
Denoting by C the complex conjugation Cip = ip, one immediately finds 
Proposition 3 The Hamiltonian (3.10) is CT and PCT invariant, i.e. 

[H, CT] = [H, PCT] = . (4.17) 

Therefore, if 

HiP e = Eip E , (4.18) 

then 

HF^pE] = EF[i^\ . (4.19) 
Clearly, CT invariance does not produce any conserved quantity since CT is anti-unitary. 
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5 Complex eigenvalues 

Interestingly, H being self-adjoint admits generalized eigenvectors with complex eigenvalues 
[19, 23, 20, 21]. Let f$ be distributions satisfying 

£/o" = 0, F/o + = 0. (5.1) 

Clearly, 

fo(x) = 6(x), g(x) = l. (5.2) 

Its easy to see that 



Hf± = ±ip±. (5.3) 



Let us define two families: 



/-:=1JL^/-, /+:= £»/+. (5.4) 

vn! vra! 



/"(*) = i-pf *<»>(*) , /+(*) = -= . (5.5) 
yn! v«' 



One finds 



Moreover, 

Hf± = ±E n f±, (5.6) 

where 

£ n :=z 7 (n+^) . (5.7) 



Clearly, both f n and / n are tempered distributions, i.e. € S' . Evidently, they are related 
by the Fourier transformation: 

jn 

F[f+] = V^i n f~ , F[f~] = —f+ . (5.8) 



Let us observe, that these two families of generalized eigenvectors have two remarkable prop- 
erties: 

/oo 
fn ( x ) fm ( x ) dx = 5 nm , (5.9) 
-oo 

and 

oo 

Y J fi{x)f-{x') = 8{x-x') . (5.10) 

n=0 



8 



These formulae remind one of the basic basic properties of proper (Hilbert space) eigenvectors: 
if A is a self-adjoint operator in TL and 

Mk = AfcV'fe , (5-11) 
where ipk are normalized vectors in Ti., then 

J 1pn(x)lp m (x) dx = 5 nm , (5.12) 

and 

^^ip n (x)ip n (x') dx = S(x — x') . (5.13) 

n 

Obviously, there is no complex conjugation in (5.9) and (5.10) since are real functions. 
Now, for any <p G Z one has 

= E ^^(-1)"^ = Y.fn (*)</» 10) • (5.14) 
On the other hand, for any <fi £ V, its Fourier transform F[(f>\ G 2, and 

= ^ / e fe fM(fc)rffc = -^=| e **^ F[ ^" )(0) (-l)"fc"dfc 
= ^F[/+](x)(/-|FM) = ^F[/+](x)(F[/-]|0) 

n n 

= E/n(*)</nl^> " (5-15) 
n 

Hence, we have two spectral decompositions: 

l^>=El/n)</nl^> ^ Z, (5.16) 

n 

and 

W) = J2\fn)(fn\4>) m (5.17) 

n 

In section 7 we derive (5.16) and (5.17) from the spectrum of H. So let us look for cr(H). 

6 Spectrum 

The Hamiltonian (3.10) has a continuous spectrum cr(H) = (— oo, oo). Since, the Hamiltonian 
(3.10) is parity invariant each generalized eigenvalue E G R is doubly degenerated: 

J?V± = E%l>g . (6.1) 
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The above equation may be rewritten as the following differential equation for 

x^(x) + \) ^ • (6-2) 

To solve (6.2) let us introduce the following distributions [13] (see also [24]): 
A / x x x>0 a/0 x>0 



x < ' ~~ ' [ |x| A x < 

with A € C (basic properties of x± are collected in the Appendix). It is, therefore, clear that 
the generalized eigenvectors V'i may be written as follows: 

^f(x) :=-J= x^ +1 /2). (6 . 4) 

It turns out that are well defined tempered distributions for all Actually, instead 
of ip± one may work with eigenvectors of the parity operator P: 

V>fven = ^ K + ^) , (6-5) 

^dd = ^ W - ^) • (6-6) 

Obviously 

P V>fven = Y^en , P ^dd = "V&d ■ (6-7) 

These distributions of definite parity are given by: 

Km = ^= \X\-W*V*> , Vo^dd = S ign(x)|x|-^/^ 1 /2) j (6 . 8) 
2v/7T7 2^/7T7 



(see [13] and [24] for the properties of |x| A and sign(x)|x| A ). 

With the normalization used in (6.4) one proves [25] orthonormality: 



Vl 1 (x)4>± (x) dx = 5{E l - E 2 ) , (6.9) 

and completeness: 



V'f (x)^±(x') dE = 8{x - x') . (6.10) 
Therefore, due to the Gelfand-Maurin spectral theorem one has 

<K*)=E / dE^(x)(^\cf>) , (6.11) 
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for any <f> € S, and the corresponding spectral resolution of the Hamiltonian has the following 
form: 



dEE\^)(^\ . (6.12) 

± J 

There is another family of energy eigenvectors directly related to ij>±. Due to (4.19) one has: 

HFty± E ]=EFty± E ] . (6.13) 
The Fourier transform of ip± is given by (cf. [13] and the Appendix): 

F[r± E ]{k) = ±^= exp [±| - £)] r + ±) (fc ± ^0)-^/^V2) . (6 . 14) 

One shows [13] that i^i] are well defined tempered distributions for any E € R. Moreover, 

y FtV-l 1 ]^) *Vf]0«0 = 5(E 1 - E 2 ) , (6.15) 

and 

J F[i> E ]{x) F[i/>£](x') dE = 5(x - x 1 ) . (6.16) 

Hence, following the Gelfand-Maurin theorem, we have further spectral decompositions: for 
any tp £ S 

V>(*) = W dEF[^ ± E ]{x){F[^ ± E ]\^) , (6.17) 
± J 

and for the Hamiltonian itself: 

• E I dEE\F[iP ± E ])(F[^ E ]\ . (6.18) 
± J 



7 Analyticity of energy eigenvectors 

Let us continue the energy eigenvectors ip E and F[if)± E ] into the energy complex plane E 6 C. 
It turns out [13] (see also the Appendix) that ip± has simple poles at E = —E n , whereas 
F[tp± E ] has simple poles at E = +E n , with E n defined in (5.7). Therefore, the poles of 
energy eigenvectors considered as functions of the complex energy correspond exactly to the 
complex eigenvalues of H which we found in Section 5. One easily computes the corresponding 
residues: 

Res(Vf (*); -E n ) = i ( T 1)-J^ S -^M , ( 7 .1) 



and 



Res(F[^ E (x)]; +E n ) = ±^ T i)^L^l x n . (7 . 2) 



11 



Hence, residues of ip± and F[ip ± E ] correspond, up to numerical factors, to the eigenvectors 
7^(5.5): 



Res(^(x);-E n ) ~ /" , (7.3) 



and 



Res(F[iP± E (x)];+E n ) ~ /+ . (7.4) 
Any function <fi G S C L 2 (R) gives rise to the following functions of energy: 

R3 E — ► (^±10) G C , 



and 



- <F[^]|0>eC. 



Let us introduce two important classes of functions [26]: a smooth function / = f(E) is in 
the Hardy class from above Tt\ (from below TL 2 _) if f{E) is a boundary value of an analytic 
function in the upper, i.e. lvuE > (lower, i.e. ImE < 0) half complex E'-plane vanishing 
faster than any power of E at the upper (lower) semi-circle \E\ — > oo. Now, define 



:= {0 e «S| <V±|0> , (7-5) 



and 



d> + :={</» G s\(Fty± E }\<f>) en 2 + } . (7.6) 
Proposition 4 <I> + n = {0}. 

Proof. Clearly, if <j) G then ( ^±10 ) is a smooth function of E G R. Suppose, that G $+, 
that is 

(F[^ E m = m E \m)^'Hi- (7.7) 

Now, due to the Paley- Wiener theorem [12] the inverse Fourier transform of F[<p] 

1 r 00 

F- l [F[4>]](E) = ^J jmt)e- itE dt , (7.8) 

vanishes for E > 0. Therefore, ^(-E) = for E > 0, and hence /(E) cannot be a smooth 
function of E. □ 
Our main result consists in the following 

Theorem 2 For any G 3>± one /ias 

0-(x) = ^/-(x)(/+|r) , (7.9) 

n 

and 

+ (x) = £/+(*)< /-|</> + > • (7.10) 
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Proof. Due to the spectral formula (6.11) one has, for cf> £ C 5: 

/oo 
dE^ix)^-) . 
-oo 



(7.11) 



Now, since {il>±\4> ) G Ti?_, we may close the integration contour along the lower semi-circle 
\E\ — ► oo. Hence, due to the residue theorem one obtains 



± n 



E—E n 



(7.12) 



Using the definition of ^ 



one finds 



1 



x^- lEh+l/2) 4>-(x)dx , (7.13) 



(x) dx . 



(7.14) 



Therefore, inserting into (7.12) the value of the residue given in (7.1) one gets finally 

0-( x ) = j-^M I [(-i)"^ + ( f ) -( x )dx = J2(- 1 ) n ^r- [x n <f>-(x)dx 

Tl. J Tl. J 

n n 

= £/n(*)(/n + l<n- (7.15) 



To prove (7.10) let us use another spectral formula (6.17): for any <p + £ & + C S 

/oo 
dEF[^ E }(x)(F[^ E }\r) . 
-oo 



(7.16) 



Now, since {F[il> ± E ]\(j) ) G 7i\, we may close the integration contour along the upper semi- 
circle \E\ — > oo. Hence the residue theorem implies 



+ 2 ™ £ £ ResOF[Y>i E (x)]; +£ n ) ( Ffo 



± n 



(7.17) 



Now, using once more the formula for ip± one finds 



<F[x^+> . 



(7.18) 
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Hence, inserting the values of residues (7.2) into (7.17) and using the formula for (see 
(A. 9)) one has 

V27T n n\ L 



= ^E(" 1 ) n - / (-i) n+1 FK](k)-i n+1 F[x!}(k) <P + (k)dk 
= Jj](-l) n — y \{-i) n+1 i n -i n+1 (-i) n 1 [ 5^(k)ci) + {k)dk 

= E^X/nl^). (7-19) 
n 

which ends the proof. □ 
This way we have recovered (5.16) and (5.17). It is not surprising, due to the following 

Proposition 5 _ = Z and <£ + = V. 

Corollary 1 We have two spectral decomposition of H : 

H = Y,E n \f-){ft\ on (7.20) 



and 



H = Y J E n\ti)(fn 



on $ + . 



(7.21) 



8 Resonances and the quantum damping 

Finally, let us turn to the evolution generated by the Hamiltonian (3.10). Obviously, it 
generates a 1-parameter unitary group 



U(t) = e 



-iHt 



on L 2 (R). It follows from (4.3) that 



ipt(x) = U(t)fp(x) = e 7t/2 V'(e 7t x) . 



(8.1) 



.2) 



The above formula is well defined for any t £ R and clearly, as we already showed, the theory 
is time-reversal invariant: if ip(t) solves the Schrodinger equation so does Tip(t) := ip(—t). 
Therefore, working in L 2 (R) we do not see any damping at all. Now, let us construct two 
natural Gelfand triplets: 



$± C L 2 (R) C $' ± 



.3) 



If G then 



-iEt I <E\ 



(8.4) 
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Hence 4>~(t) G <£„ only for t > 0. Similarly, if <p + G <I> + , then <^ + (t) G <J?+ only for t < 0. 
Therefore, the restriction of the unitary group ^7(t) on L 2 {M) to <I>± no longer defines a group. 
It gives rise to two semigroups: 

U-(t) : — ► , for t>0, (8.5) 

and 

U+(t) : <$>+ — > $+ , for t < . (8.6) 
Due to (7.20) and (7.21) one has: 

fW=^)f = E e_7(n+1/2) *l/n)^nlr) , (8-7) 
n 

for i > 0, and 

= U(t)4> + =Y J ^ n+mt \fn){fn\<l> + ) , (8-8) 

for i < 0. We stress that 07" does belong to L 2 (R) also for t < (i > 0). However, 

07 G G $+) only for t > (t < 0). This way the irreversibility enters on a purely 

Hamiltonian level by restricting dynamics to the dense subspace §± of I? (R) . 

Clearly, formulae (8.7) and (8.8) are quantum analogues of the classical damping laws: 

x(t) = e'^x , t > , (8.9) 

and 

p{t) = e + ^p , t < . (8.10) 

Finally, let us recall that the time reversal operator T establishes an isomorphism between 
<5_ and <E> + . Therefore, each solution 

^ = U.(t)tf>- , (8.11) 

with <J)~ G 3>_ is mapped into 

T(^) = U-(-t)T(</>-) = f/+(W<n , * < . (8.12) 
Conversely, any solution 

<j>t = U + {t)ct> + , (8.13) 

with G 3>+ is mapped into 

T(0+) = l/ + (-t)T(^+) = U-(t)T(<f>+) , t > . (8.14) 

Summarizing, quantum dynamics is irreversible on <!>_ and <I> + . This irreversibility is caused 
by quantum damping, or, equivalently, by the presence of resonant states (5.5). It should 
be stressed that it is not an energy that is dissipated. Clearly, the Hamiltonian (3.10) can 
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not be interpreted as a system energy — it was used to define a Hamiltonian dynamics of the 
enlarged system on L 2 (IR). The quantum damped system is not defined on the entire L 2 (R) 
but rather on a dense subset 3>_ C L 2 (R). As we saw it imposes the restriction upon the 
time evolution such that it is defined only for positive t. A quantum damping may be seen 
as follows: let 6q £ be an initial state then the probability density for a particle position 
evolves in time as follows: 

pt(x) = \6^(x)\ 2 = e^o (e 7t *)| 2 = e^o(e^x) , (8.15) 
and hence in the limit t — ► +00, one finds pt(x) — ► S(x). Indeed, for any e > 

/e ree jt 
Pt(x)dx = / po{x)dx — ► 1 , (8.16) 
-e J -eel 1 

for t — ► +00. Clearly, it corresponds to the classical behavior x(t) = e~ yt XQ — > 0. 

In a forthcoming paper we are going to show that also more complicated damped systems, 
e.g. the damped harmonic oscillator, give rise to irreversible dynamics. 



A Appendix 

The regular tempered distribution x\_ (with A G C) given by 



POO 

6\x\) := \ x x <p{x)dx , 
Jo 



(A.l) 



for any <fr G S, is well defined for Re A > —1. However, it may be easily extended to the 
region ReA > — 2 due to the following regularization formula: 



poo pi roc 6(0) 

/ x x 6(x)dx = / x x [6(x) - 6(0)]dx + / x x 6(x)dx + s 
Jo Jo Jl 



A + l 



(A.2) 



which holds for A / —1. In the same way one may extend the distribution x+ to the region 
Re A > — n — 1 using the formula 



x x 6(x)dx 



(x) - 6(0) - x6'(0) - 



x 



n— 1 



(n-1)! 



(n-l) 



(0) 



dx 



+ 



/ x x 6(x)dx + 
J 1 7 — 1 



^"^(O) 



(A.3) 



which holds for A / —1, —2, . . . , — n. The above formula shows that ( 6\x\ ) as a function of 
A G C has simple poles at A = —1, —2, . . ., and the corresponding residue at A = — k equals 
6^-^(0) /(k-\)\. 
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Using the same arguments one shows that the distribution x A may be extended to the 
region Re A > — n — 1 via: 

/0 poo 
x x <fi(x)dx = / x X (j){—x)dx 
-oo JO 

roc r (_-\)n—l T n—l 
= / X X 0(- X )-0(O)+X0'(O)-...-^ ^L_^("-1)(0) 

A L (ra — 1J! 

which holds for A 7^ —1, —2, . . . , — n. Hence, ( 4>\x^_ ) has simple poles at A = —1, —2, . . ., and 
the corresponding residue at A = —k equals (—l) k ~ 1 (/)( k ~ 1 \o)/(k — 1)!. 
The Fourier transforms of x± 

F[xi](k) = ^= [ e ikx x\dx , (A.5) 

V27T J 

are given by the following formula [13] 

F[xi]{k) = ±-Le ±iA "/ 2 r(A + l)(k + iO)-^ 1 , (A.6) 

V27T 



where (k ± i0) a is a distribution defined by: 

(fc ± i0) a = + e ±ia7r ^ . (A.7) 

Due to the Euler T-function the formula (A.6) has single poles at A = —1, —2, . . .. Note, that 
although both /c° and k* have poles at a = —1, —2, . . ., the distribution (k ± i0) a is well 
defined for all a G C. Indeed 

lim (jfe ± iO) Q = lim (fe+ + (-l) n fc°) , (A.8) 

and, due to (A. 3) and (A. 4), the singular parts of k" and fc°, at a = — n, cancel out. In 
particular, for A = n G N, one obtains (cf. [13]) 

F[a£] (fc) = -L [(±i) n+1 n! fc-™- 1 + (Ti) n 7r<f < n > (jfe)l . (A.9) 

V 27T L 
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